ABSTRACT
A thorough understanding of the inhalation dynamics of infectious aerosols indoors and infection dynamics within the host by inhaled viruses such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) plays an important role in the assessment and control of infection risks indoors. Here, by combining computational fluid–particle dynamics (CFPD) and host–cell dynamics (HCD), SARS-CoV-2 infection dynamics in the mucus layer of the human upper airway were studied. To reproduce the diffusive and convective transport of the virus in the nasal cavity–nasopharynx by mucociliary motion, a three-dimensional (3D)-shell model with a mucus layer was developed. The initial virus concentrations for HCD calculation were estimated based on the deposition distribution of droplets with representative sizes analyzed by CFPD. To develop a new HCD model, the target-cell-limited model was integrated with the convection–diffusion equation. Additionally, the sensitivity of the infection rate β to the infection dynamics was systematically investigated. The results showed that the time series of SARS-CoV-2 concentration in the mucus layer strongly depended on diffusion, convection, and β. Although the SARS-CoV-2 dynamics obtained here have not been verified by corresponding clinical data, they can preliminarily reveal its transmission mode in the upper airway, which will contribute to the prevention and treatment of coronavirus disease 2019.
ABSTRACT
Multilayer diffusion problems have found significant importance that they arise in many medical, environmental, and industrial applications of heat and mass transfer. In this article, we study the solvability of a one-dimensional nonhomogeneous multilayer diffusion problem. A new generalized Laplace-type integral transform is used, namely, the Mρ,m-transform. First, we reduce the nonhomogeneous multilayer diffusion problem into a sequence of one-layer diffusion problems including time-varying given functions, followed by solving a general nonhomogeneous one-layer diffusion problem via the Mρ,m-transform. Hence, by means of general interface conditions, a renewal equations’ system is determined. Finally, the Mρ,m-transform and its analytic inverse are used to obtain an explicit solution to the renewal equations’ system. Our results are of general attractiveness and comprise a number of previous works as special cases.